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Abstract 

The electric dipole form factor of the nucleon stemming from the QCD 9 term 
is calculated in chiral perturbation theory in leading order. To this order, the form 
factor originates from the pion cloud. Its momentum dependence is proportional 
to a non-derivative time-reversal-violating pion-nucleon coupling, and the scale for 
momentum variation — appearing, in particular, in the radius of the form factor — 
is the pion mass. 
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Electric dipole moments (EDMs) of the neutron and of atoms have long been of interest 
because precise experiments are possible, providing strict bounds on possible sources of 
CP violation. (For a review of both experimental and theoretical results, see for example 
Ref. [1].) New experiments with ultracold neutrons in preparation at various laboratories 
(LANL, PSI, ILL, Miinchen) promise to provide even higher precision for the neutron 
EDM than the current value dn = (—1.0 ±3.6) ■ 10^^^ e cm [2]. While neutron experiments 
yield bounds on the neutron EDM directly, atomic effects are sometimes sensitive to the 
nuclear Schiff moment, which in turn receives a contribution from the radius of the nucleon 
electric dipole form factor (EDFF). 

Nucleon properties are notoriously difficult to calculate directly from QCD. At low 
momenta, Q ~ -C Mqcd, where 171^^ is the pion mass and Mqcd ~ 1 GeV is the 
typical mass scale in QCD, they can nevertheless be described in terms of an effective 
field theory (EFT) involving nucleons and pions (and delta isobars), known as chiral 
perturbation theory (^PT). (For an introduction, see for example Refs. [3, 4].) In this 
theory, long-range effects due to the light pions are separated from short-range effects 
due to all higher-energy degrees of freedom. Observables are systematically expanded 
in powers of Q/Mqcd (times functions of Q/m,,). The pion cloud around the nucleon 
can generate contributions that are non-analytic in momenta and quark masses, and 
are calculable in terms of pion and nucleon (and delta) parameters. In contrast, short- 
distance physics contributes powers of and that can only be calculated using 
detailed knowledge of the QCD dynamics. 

The leading non-analytic contribution to the nucleon EDM coming from the QCD 9 
angle was first calculated in xPT in Ref. [5], and rederived many times since, see e.g. 
Ref. [6]. A measurement of the neutron EDM cannot disentangle this long-range effect 
from short-range, analytic terms. Since a cancellation between different functions of 
is unlikely, the leading non-analytic contribution serves as an estimate of the EDM. 

Here we consider the full nucleon EDFF in leading order in xPT. The other nucleon 
form factors have all been calculated previously in this theory at the first few orders: 
electric [7], magnetic [7], and anapole [8]. Needless to say, the momentum dependence of 
the EDFF is not as easily accessible experimentally However, it might become of interest 
after the neutron EDM (the EDFF at = 0) is measured. As we will show below, the 
EDFF is in leading order a calculable isovector function of Q^/m^, entirely determined 
by the pion cloud. It is that determines the scale of variation in the form factor — an 
example of a low-energy theorem. In particular, the radius of the EDFF is fixed by m^, 
as first derived in Ref. [9]. If measured, the neutron EDFF would allow a determination 
of the T- violating pion-nucleon coupling (within the precision of the xPT expansion). 

At Q ~ m^, the nucleon is essentially non-relativistic, because mjv ~ Mqcd- Pions 
must explicitly be accounted for in the theory, since they are the fight (pseudo-)Goldstone 
bosons corresponding to the spontaneous breaking of (explicitly-broken) chiral SUl{2) x 
SUr{2) ~ S0{4) symmetry down to SUl+r{'2) ~ 5*0(3). Also, since the mass difference 
between the A isobar and the nucleon is comparable to the pion mass, m/\ — trn ^ 
2m^, the A should be included as well (which is not much of a problem since at these 
momenta, the A is also a nonrelativistic object). This can be done by first constructing 
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the most general Lagrangian involving nucleons, pions, and deltas that transforms under 
the symmetries of QCD in the same way as QCD itself. Along with this, one needs a 
power-counting scheme so that interactions can be ordered according to the expected size 
of their contributions. The strong-interaction Lagrangian contains an infinite number of 
terms that can be grouped according to the index A: 

oo 

, A^d + //2-2, (1) 

A=0 

where d is the number of derivatives, powers of m„ and/or powers of rriA — ttin, and / is 
the number of fermion fields. Electromagnetic interactions are proportional to the small 
charge e, and it is convenient to account for factors of e by enlarging the definition of d 
accordingly. 

The technology for constructing the Lagrangian is well known, see for example Ref. 
[3] . We will need the following T-conserving terms, which either obey chiral symmetry or 
break it in the same way as the quark-mass terms do: 

^%em = Id.tt . D^-K - ^TuW + Niv ■ DN - ^N(S ■DT.7r)N + ... (2) 

Here tt denotes the pion field in a stereographic projection of SO{A)/ S0{3), with D = 
1 -|-7r^/4/^ and f^^ = 93 MeV the pion decay constant, N = [p n)^ is a heavy-nucleon field 
of velocity f'* and spin 5^ [S^ = (0,(7/2) in the nucleon rest frame where v*^ = (1,0)), 
Dn is the covariant derivative, with {D^^^)ab — D~^{5abd^ + iee^ab^n) for a pion and 
D^/^) = d^ + iT ■ {tt X D^-k) /Afl — ieAn{l + t^) /2 for a nucleon, and ". . ." stands for other 
interactions with more pions, nucleons and/or delta fields that are not explicitly needed 
in the following. Note that the pion mass term comes from the explicit breaking of chiral 
symmetry by the average quark mass fn = (m„ + md)/2, so = 0{Mqcd^)- Isospin 
breaking from the quark mass difference, which is proportional to me — 0{m'^e/MQCD) 
with e = {md — mu)/{mu + md) — 1/3, does not appear at this order. Note also that 
at this order the nucleon is static and couples only to longitudinal photons. Kinetic 
corrections and magnetic couplings have relative size 0{Q / Mqcd) and appear in /^s2/em- 
The same is true for the delta isobar, including the nucleon-dclta transition through 
coupling to a transverse photon. The pion-nucleon coupling in Eq. (2) is not determined 
from symmetry, but is expected to be C(l), and indeed qa — 1.267. The Goldberger- 
Treiman relation gA'^N = fngnNN holds in lowest order. A term in C J^i^^ provides an 
Oiim-KjMQCD)^^ correction that removes the so-called Goldberger-Treiman discrepancy. 

In addition, we need P- and T- violating interactions. The exact form of these inter- 
actions will depend on the mechanism of CP violation. Just above Mqcdi T-violating 
interactions among quarks and gluons can be classified through their scale dimensions, 
starting with the 6 term. In a basis where the quark fields q — {u d)^ have been appro- 
priately rotated [10], we can write 

C^^^ = mj qi'f^q + (3) 
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where = mumd/ {rriu + = fh{l — £:^/4)/2 ~ fh/2. Here ". . ." represents higher- 
dimension operators — such as the quark EDM and color EDM, the Weinberg three- gluon 
operator, and four-quark interactions — which wc will neglect. 

The 9 term generates T-violating interactions in the low-energy EFT. It is the fourth 
component of the 5*0(4) vector P = {qTq,qi'j5q), so it generates EFT interactions that 
transform as T-violating, fourth components of 50(4) vectors made out of hadronic fields. 
The lowest-dimension operator of this type is 

£W = -gATr-TTA^ (4) 

where (^o is an / = T-violating ttNN coupling. From dimensional analysis, go ~ 
0{mj)/ f-jj) = 0{ml9 / f-^MqcD)- Because of the factor of implicit in ^o, i^f^ ef- 
fectively has index A = 1. 

Interactions with higher indices can be constructed in similar ways. T-violating nucleon- 
dclta transitions through pion emission, for example, involve at least one derivative, and 
are thus suppressed by at least 0{Q / Mqcd) relative to those stemming from the La- 
grangian (4). Among the higher-order operators, particularly relevant here are short-range 
contributions to the nucleon EDM, 

= N{do + d^r^){S^v, - S,v^)NF^\ (5) 

where do {di) is a short-range contribution to the isoscalar (isovector) EDM of the nucleon. 
Prom dimensional analysis, — 0{emifi / Mq^j-,) — 0{em'^9 / Mqqj^) . Direct short-range 

contributions to the momentum dependence of the EDFF first appear in being 
further suppressed by C((Q/Mqot)^)- 

Here we denote by iJ'^^ the T-violating nucleon current that interacts with the elec- 
tron current —ieej^e via the photon propagator iDj^^, = —i{r]^y/q^ + . . .) to produce a 
contribution 

iT = -iee{k')^>'e{k) D,,{q) N{p')rM)N{p), (6) 
to the electron-nucleon S matrix. We have q^ = {p — p'Y = —Q^ < 0. We can write 

JM - 2 {F^^\-q') + F^j^\-q')T,) {S^v -q-S- qv^), (7) 

where F^\q'^) {Fj^\Q^)) is the isoscalar (isovector) EDFF of the nucleon, with Fj^\o) = 
di the corresponding EDM. 

When we consider the sizes of specific contributions to the neutron EDFF, the first 
tree-level contribution comes from the vertex generated by the Lagrangian of Eq. (5). The 
lowest-order one-loop graphs are built out of one vertex from the T-violating Lagrangian 
(4) and all other vertices from the strong Lagrangian (2). The one-loop diagrams that 
a priori could contribute to the EDFF are shown in Fig. 1. These one-loop diagrams 
contribute 0{e%fT^Q/ {Anf^^y) to Jedio): ^-^d thus give rise to an EDFF of the same order 
as the short-range, tree-level contribution to the EDM, 0{em-ifi / Mqqj^) . 
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(a) (b) (c) 



Figure 1: One-loop diagrams contributing to the nucleon electric dipole form factor in 
leading order. Solid, dashed and wavy lines represent nucleon, pion and photon, respec- 
tively; circles and squares stand for interactions from ^^g^r/em ^^'^ \ respectively. For 
simplicity only one of two possible orderings is shown here. 



Because of the isospin structure, the contribution from Fig. 1(c) vanishes. We evaluate 
diagrams in Fig. 1 at v-q = 0, as a consequence of the static nature of the nucleon in lead- 
ing order. Diagram 1(b) then also vanishes. Diagram 1(a) gives a non-zero contribution 
only to the isovector component of the nucleon EDFF. 

The contributions from Fig. 1 contain divergent pieces that can be absorbed by redefin- 
ing the tree-level contribution to the EDFF. Taking the tree and one-loop contributions, 
we have 

do = do (8) 

/i 



di — di + 



A + 2 In ■ 



(9) 



where /i is the renormalization scale introduced by dimensional regularization, and A = 
2/s — 7e + ln47r {e — 4 — d, with d the spacetime dimension). Using the Goldberger- 
Treiman relation and setting // to mjv, the piece in Eq. (9) that is non-analytic in can 
be written as {egj^NNgo/ 4:7i'^mN)\n{mN / 171^,) , in agreement with the result of Crewthcr et 
al. [5]. The short- and long-range contributions are in general of the same size, but a 
cancellation is unlikely due to the In m,r dependence of the pion contribution. Note that 
short- and long-range physics cannot be separated with a measurement of the neutron 
EDM dn^ do - di alone. 

Because the form factor is given by lowest-order loop graphs, it depends on the com- 
bination Q'^/m\ only. Because the only non- vanishing contribution to this order comes 
from Fig. 1(a), the combination is actually (5^/(2m^)^. The EDFF defined in Eq. (7) is 
found to be given by 

F^n\Q') = do (10) 
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Figure 2: The function F((5^/(2m7r)^) that enters the isovector electric dipole form factor 
F^^^ xPT in leading order, Eq. (12) (solid line); and quadratic approximation, Eq. (13) 
(dotted line). 

where 



n.).3^.,/l..ln|.Iii|±i)-i;. (12) 

Note that F{Q) = 0, so indeed this does meet the standard of being a form factor. The 
function F is a testable prediction of xPT: it is plotted as a function of in Fig. 2. 
For X <^ 1 we can expand F{x) in powers of x, 

F{x) =x + 0{x'^). (13) 

This approximation is compared to the full functional form (12) in Fig. 2. The variation 
of the form factor with Q can be characterized at very small momenta by the isovector 
EDM square radius (defined in analogy to the charge square radius). 

While the EDM vanishes in the chiral limit, the radius of the EDFF is finite. As was 
noted in Ref. [9], S' = {rld)^^^^ = ^QaQo/ ^'&t^'^ f^nin'^ is the leading electromagnetic 

contribution to the Schiff moment of the the proton (neutron). 
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The mass scale governing the variation of F{Q'^/{2m'^)) is in fact 2m^, since it is 
generated by the pion cloud. The relative importance of the variation of F^\q^) 
depends in addition on the ratio egAgo/^'^T^^ fi^di. This ratio is a constant in the chiral 
limit. Since we expect it to be C(l), it is likely that m.^ sets the scale of the variation 
of Fj^\Q'^) as well. Unlike other nucleon form factors, a dipole approximation with a 
mass scale close to the rho-meson mass, m^, should not be good even as a numerical 
approximation for the EDFF. 

All other contributions to the EDFF are formally of higher order. Next-order effects 
come from diagrams with one insertion of an operator with one more derivative. First, 
there are operators in 



+1 (1 + «o + (1 + «i)t3) vI'S^NEP'^ + . . j + . . . (15) 

These are nucleon recoil corrections to the static limit and photon coupling to the nucleon 
magnetic moment (including the isoscalar and isovector anomalous magnetic moments, 
kq — —0.12 and ki — 3.7 respectively). Second, there is a new operator in 

= • D^TZ NS^'N, (16) 

with an undetermined coefficient ho = 0{mi,6 / f'^Mqco) = ^ij^V^ / fl^QCo)- The corre- 
sponding contributions to the nucleon EDFF are also calculable, as they are still given by 
long-range physics associated with the pion cloud. Although the strong-interaction cor- 
rections depend on known parameters, Eq. (16) introduces a new T-violating parameter. 
Both types of contributions — from Eqs. (15) and (16) — are suppressed by 0{Q / Mqcd) 
relative to the contributions calculated here. 

Note that if £^ contained other non-derivative, T-violating pion-nucleon interactions, 
TTsNN and tt^NtsN [11], our previous results would remain unchanged. However, the 
chiral-symmetry transformation properties of the 9 term do not yield such terms. 

In conclusion, we have presented results for the electric dipole form factor of the nucleon 
due to the 6 term in leading order in xPT. We have shown that the the EDFF is isovector, 
with a dependence determined by a non-derivative T-violating pion-nucleon coupling 
and the pion mass, see Eqs. (11) and (12), and Fig. 2. Under the assumption that 
higher-order results are not afflicted by anomalously-large dimensionless factors, the error 
of our results at momentum Q should be ~ Q/rUp. 
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